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ABSTRACT 

We present three years of long-cadence and over one year of short-cadence photometry 
of the central star of the Planetary Nebula NGC 6826 obtained with the Kepler space- 
craft, and temporally coinciding optical spectroscopy. The light curves are dominated 
by incoherent variability on time scales of several hours, but contain a lower-amplitude 
periodicity of 1.23799 d. The temporal amplitude and shape changes of this signal are 
best explicable with a rotational modulation, and are not consistent with a binary 
interpretation. We argue that we do not observe stellar pulsations within the limi- 
tations of our data, and show that a binary central star with an orbital period less 
than seven days could only have escaped our detection in the case of low orbital in- 
clination. Combining the photometric and spectroscopic evidence, we reason that the 
hourly variations are due to a variable stellar wind, and are global in nature. The 
physical cause of the wind variability of NGC 6826 and other ZZ Leporis stars is likely 
related to the mechanism responsible for wind variations in massive hot stars. 

Key words: stars: variables: other - stars: early-type - stars: mass loss - stars: 
rotation ~ stars: winds, outflows - planetary nebulae: individual: NGC 6826 



1 INTRODUCTION 

The ZZ Leporis stars (Handler 2003) are a group of fourteen 
variable Central Stars of Planetary Nebulae (CSPN), named 
after their prototype, the central star of IC 418. The cause of 
their variability remains to be fully understood. Photometric 
variations are present on two time scales, of the order of days 
and of the order of hours. The most plausible mechanisms 
for the variability are pulsation or variations in the stellar 
mass loss rate - or both. Binary-induced light variations can 
in most cases be excluded due to the observed hght curve 
shapes, the lack of a dominant periodicity, and the lack of 
corresponding radial velocity changes. 

The reason why pulsations may be expected to be 
present in ZZ Lep stars is evident from their positions in 



the HR diagram. They are located at the intersection of 
the post- Asymptotic Giant Branch evolutionary track with 
the instability strip of the fi Cephei pulsators. Gautschy 
(1993) showed that stellar models reminiscent of those stars 
have pulsationally unstable eigenmodes. Zalewski (1993) 
confirmed these results and showed that nonlinear mod- 
els in the temperature range of the ZZ Lep stars undergo 
quite complicated hght and radial velocity variability on 
time scales of several hours. 

Wind variability, on the other hand, was proposed as 
a mechanism by Mendez, Verga & Kriner (1983). UV spec- 
troscopic studies (e.g., Prinja et al. 2012) support this idea, 
because all ZZ Lep stars examined indeed showed clear ev- 
idence for variable winds. It is believed that the temporal 
behaviour of the winds of those CSPN is similar to that of 
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massive OB stars (e.g., Kaper et al. 1997), is related to the 
stellar rotational period, and may be rooted close or at the 
stellar surface (e.g., FuUerton et al. 1997). Further strong 
support for this interpretation was recently presented by 
Prinja, Massa & Cantiello (2012), who also argued that the 
wind variations of the central star of NGC 6543 could be 
causally connected to sub-surface convection. The radiation 
pressure driven winds of hot stars generally appear to exhibit 
extensive variability in their stellar winds, whether on the 
main sequence, as supergiants, Wolf-Rayet stars or CSPN. 

Whatever the cause of variability of the ZZ Lep stars, 
it will have implications for our astrophysical knowledge of 
the evolution of CSPN. They cross the HR diagram on time 
scales of a few thousand years under the influence of mass 
loss. If ZZ Lep stars pulsate, measurements of oscillation 
period changes could closely constrain their evolutionary 
speed. Model calculations (e.g.. Blocker 1995) imply that 
this speed is highly mass dependent and therefore CSPN 
masses could be accurately determined. If the complicated 
light variations of ZZ Lep stars were due to multimode pul- 
sations, they might even be accessible to asteroseismology. 
Constraints on CSPN mass loss behaviour, on the other 
hand, are also an asset for understanding their evolution, 
spectral characteristics and even for the shaping of the sur- 
rounding nebulae (e.g., Huarte-Espinosa et al. 2011). 

One of the reasons why the variability mechanism of the 
ZZ Lep stars is so hard to pinpoint may be their complicated 
behaviour in combination with a lack of supporting data. 
Ground-based observations suffer from gaps due to daylight 
interruptions (even in case of multisite campaigns), whereas 
time on rather large telescopes or UV spacecraft required 
for spectroscopic monitoring is very limited. In addition, ZZ 
Lep stars are rather faint and surrounded by bright nebulae 
that complicate measurements and their interpretation. On 
the other hand, their variations occur on shorter time scales 
than those of massive OB stars, which are of the order of 
one day, which can be an observational advantage. 

One ZZ Lep star is located within the field of view of 
the Kepler mission (Koch et al. 2010). HD 186924 = KIC 
12071221 {V = 10.4) is the central star of NGC 68243- Its 
photometric variability was first reported by Bond & Ciar- 
duUo (1990) and confirmed by Handler (1998) and Michalska 
(2000). The light range of the variations was several hun- 
dredths of a magnitude, and the time scale several hours, 
with no clear evidence for a periodicity. UV spectroscopy 
(Prinja et al. 2012) revealed the presence of recurring Dis- 
crete Absorption Components (DACs) in the absorption 
troughs of wind-sensitive lines, clear evidence for variable 
wind structure. 

Jevtic et al. (2012) examined the first ~300d of Ke- 
pler photometry of NGC 6826 by means of nonlinear time 
series analysis and phase-space reconstruction. They noted 
that the variability is coherent on time scales of up to about 
1.5 hr, but becomes consistent with noise on time scales ex- 
ceeding 10 hr. In the present work, we analyse all presently 
available Kepler time-series photometry with more tradi- 
tional techniques. Combining the photometric results with 



^ In the remainder of this paper, we will use the name of the 
Planetary Nebula also as the designation of the central star. 



constraints from optical time-series spectroscopy, we discuss 
interpretations of the central star variability. 



2 OBSERVATIONS 
2.1 Photometry 

Kepler photometry is provided in two modes, short cadence 
(SC) and long cadence (LC). In SC mode (Gilliland et al. 
2010a), a data point is obtained every 58.8 seconds, whereas 
LC mode (Jenkins et al. 2010b) yields a data point every 29.4 
minutes. The data are arranged in quarters of three months 
each, with the exception of the first two. QO consists of 10 
days of commissioning data and Ql of one month of science 
data. In between the quarters there are gaps of about one 
day for data download and spacecraft rotation. 

NGC 6826 was and still is observed in LC mode from Ql 
(starting on 13 May 2009) onwards. After examination of the 
first LC data, SC measurements were requested within the 
Kepler Asteroseismic Science Consortium (KASC, Gilliland 
et al. 2010b), which were obtained from Q7 - Qll. Kepler 
photometry is available as "corrected" and "uncorrected" 
data. In this sense, "uncorrected" means that the data have 
undergone standard reductions such as bias correction, fiat- 
fielding and correction for the smear induced by the readout, 
whereas it has been tried to remove instrumental systemat- 
ics in the "corrected" data. We started from the "uncor- 
rected" data because a comparison showed that these will 
result in a lower noise level. Both the LC and SC light curves 
were first converted to magnitudes, visually inspected, and 
in case of the presence of simultaneous LC and SC data com- 
pared. Long-term drifts were assumed to be instrumental 
(for instance, arising from differential velocity aberration) 
and were filtered out with polynomials up to third order to 
account for the instrumental drifts properly whilst avoid- 
ing to "overfit" the data in one-month data chunks. Hence, 
our analysis will not be sensitive to periodicities longer than 
about a week. Some very few outliers and bad sections of 
data, e.g. after pointing losses of the spacecraft whose times 
are documented were removed by hand. The final data set 
on which this paper is based comprises the Ql - Q12 LC 
data, which span 1051 days with a duty cycle of 93.3%, and 
the SC data containing 468 days of observation with a 92.3% 
duty cycle. We estimate a precision per data point of 0.17 
mmag for the SC data; the apparent precision of the LC data 
however only is 0.43 mmag per point, implying that most of 
the "noise" in these measurements is actually signal. 

One difficulty with CSPN photometry is the presence of 
the nebula (that is not expected to vary on the time scales 
investigated here). Its infiuence can be suppressed by choos- 
ing filters at low nebular continuum fiux, and avoiding strong 
nebular emission lines. Another method is to mathematically 
subtract the nebular contribution or image subtraction. All 
these are not suitable for the Kepler data. We use the normal 
Simple Aperture Photometry (SAP, Jenkins et al. 2010a) 
data. Those contain some nebular fiux but are more robust 
against artificial amplitude changes coming from e.g. (quar- 
terly) pointing changes of the spacecraft than single-pixel 
data. As a consequence, the amplitudes of the light varia- 
tions reported here are smaller than they intrinsically are. 
A comparison between SAP and pixel data implies an am- 
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Figure 1. A nine-day portion of the Kepler photometry of NGC 6826. 



plitude reduction to less than one quarter of the intrinsic 
value. A subset of the SC data is displayed in Fig. 1. 

2.2 Spectroscopy 

We obtained high-resolution spectroscopy co-temporal with 
a section of the Kepler photometry. The first data set stems 
from the 3.6-m Canada France Hawaii Telescope (CFHT) 
and the ESPaDOnS spectrograph (Donati 2003). Observa- 
tions were taken on three consecutive nights from 30/31 May 
to 1/2 June, 2010 (UT), totalling 14 spectra with individual 
exposure times of 1800 s, a wavelength coverage of 3700 to 
10480 A and a spectral resolution of 68000 in "star -|- sky" 
mode. On the 2.6-m Nordic Optical Telescope (NOT), the 
FIES (Frandsen & Lindberg 1998) spectrograph was used on 
the nights of 6/7 to 7/8 June, 2010 (UT). Thirteen spectra 
were gathered with individual exposure times of 1800 s, in a 
wavelength range of 3640 to 7360 A with R = 46000. 

The optical spectrum of NGC 6826 is characterised by 
variable lines including Hell 4542 A and Hell 5411 A in 
absorption. NV 4605 A exhibits a weak P Cygni profile, and 
Hell 4686 A as well as CIV 580lA and 5812 A are in strong 
(stellar) emission. Further visible are the Balmer absorption 
wings (the line cores arc contaminated by nebular emission) , 
and NIV 638lA is a rare, but very weak, photospheric line. 
Figure 2 shows examples of the hourly variability present 
in Hell 4686 A and Hell 5411 A in the CFHT and NOT 
spectra (see Sect. 3.2). 



3 ANALYSIS 

3.1 Photometric frequencies 

The light variations of NGC 6826 are quite complicated, 
predominantly occur on time scales of a few hours and show 
no obvious periodicity (Fig. 1). We analysed the whole SC 
data set using the Period04 software (Lenz & Breger 2005). 
This package applies single-frequency power spectrum anal- 
ysis and simultaneous multi-frequency sine-wave fitting. It 
also includes advanced options such as the calculation of 



optimal light-curve fits for multipcriodic signals including 
harmonic, combination, and equally spaced frequencies. 

The amplitude spectrum of the SC data is shown in 
Fig. 3. It is dominated by a single frequency, surrounded by 
apparent noise dropping continuously in amplitude towards 
higher frequencies. No coherent signal is present between 20 
/xHz up to the Nyquist frequency of 8500 fjMz. 

However, the predominant light variations of NGC 6826 
seem to occur on this time scale (Fig. 1). Therefore, we com- 
puted Time-Fourier spectra of all the SC data. One example 
of those is shown in Fig. 4, where we chose a five-day slid- 
ing window and a time step of 0.5 days. Again, there is no 
evidence for coherent or recurrent signals. If there is power 
in the Time-Fourier spectrum, it is short-lived, as experi- 
menting with different temporal window sizes prove. In Fig. 

4 it is also visible that the amplitudes and frequencies of the 
short-term variations are not stable in time, but are variable 
as well. As there is no more astrophysical information to be 
extracted from the Time- Fourier analysis, we decided to con- 
tinue our investigations with the simplest possible methods. 

As an additional test, we divided the data set into two 
halves and computed their amplitude spectra. Their appear- 
ance is consistent with that expected from noise: their back- 
ground levels are similar, but compared to the full data set 
by a factor of about \/2 higher. This result is fully consistent 
with the earlier study by Jevtic et al. (2012), who showed 
that the variations are coherent over short time scales only. 

Turning to the apparently coherent signal, we can make 
use of the LC data and their longer time base. The ampli- 
tude spectrum of these measurements is examined in Fig. 

5 (upper panel). As in Fig. 3, there is a dominant signal, 
but its amplitude is only about 60% of that in the SC data. 
Prewhitening this signal (lower panel of Fig. 5) does not re- 
move all the power around this frequency; several peaks in 
a narrow interval (~ 0.06 /iHz) remain. In addition, a sub- 
harmonic of this signal is detected. Investigating the sub- 
harmonic with the same methodology shows a similar pic- 
ture: it is not reproducible with a single coherent frequency, 
and residual peaks in an ~ 0.09 /uHz interval remain after 
prewhitening. 
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Figure 3. Amplitude spectra of tiie siiort-cadence Kepler ptio- 
tometry of NGC 6826. 
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Figure 2. Dynamic spectra of the HcII 4686 A and Hell 5411 
A lines. The individual spectra have been ratiocd by the mean 
spectrum for each telescope. The top two panels are NOT data, 
whereas the bottom two represent the CFHT spectra. The panel 
at the very bottom shows examples of individual spectra from the 
CFHT data to illustrate the maximum extent of the line intensity 
changes. There is some correspondence between the variability 
patterns seen in these two Hell lines that is also seen in some other 
spectral features. HeH 4686 is variable to at least ±200km/s. 



To assess the significance of the detection of the two sig- 
nals, wc adopted the widely used and reliable criterion by 
Breger et al. (1993). According to this criterion, any peak 
that exceeds a signal-to-noise (S/N) ratio of four in the am- 
plitude spectrum is statistically significant. To compute the 
noise level in the presence of red noise and time-variable 
signals, wc integrated over the amplitude spectrum between 
— 8/iHz and 20 — 28/liHz and interpolated the results to the 
frequencies of interest. The dominant peak has S/N of 22.5, 
and its subharmonic of 4.4. The next tallest peak in Fig. 5 
has S/N = 3.6 and is therefore not significant. 

These findings were traced to a change of the Ught curve 
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Figure 4. Time-Fourier spectrum of one month of short-cadence 
Kepler photometry (Q9.1) of NGC 6826, with a sliding window 
size of 5 d. To accentuate the variations on short time scales, the 
dominating low-frequency signal is not shown. 



amplitude and shape over monthly time scales. Phase dia- 
grams of 2-month stretches of data, folded among the sub- 
harmonic frequency are displayed in Fig. 6. Averaging over 
such time intervals was found most suitable to suppress the 
incoherent short-term variations whilst sampling the tem- 
poral change of the periodic component of the hght curve 
well and not affecting any scientific conclusions drawn. The 
light curve evolved from a double-wave shape to a single- 
wave variation, then almost dropped to zero amplitude and 
reached almost a 0.01 mag hght range before dropping to 
an asymmetric double-wave variation again. The phases of 
light maximum and minimum varied only slightly over the 
three years of Kepler observations. 

We adopted the subharmonic frequency 9.3491 ± 
0.0001 //Hz, corresponding to a period of 1.23799±0.00001 d, 
as the fundamental time scale of the regular variability of 
NGC 6826. The reason for this choice is that in almost all 
variable stars in whose amplitude spectra a subharmonic is 
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Figure 5. Amplitude spectra of the ~ 3-year LC iCepier photom- 
etry in the low-frequency domain. The dominant signal /i does 
not cleanly prewhiten. In the lower panel, the subharmonic of /i 
also stands out. 



detected, it is the subharmonic frequency that is related to 
the physical cause of the variability. Prominent examples 
are eclipsing binaries or ellipsoidal and rotational variables. 
The only exception from this rule is the period doubling 
phenomenon, a rare and complex resonance effect in some 
pulsating stars (e.g., see Szabo et al. 2010). Aside from that, 
adopting the tallest peak in the amplitude spectrum as the 
physical variability time scale would leave the subharmonic 
unexplained. 



3.2 Spectroscopic constraints 

The 1.23799-day Kepler photometric period is not detected 
in our spectroscopic data set. However, we caution that the 
temporal sampling of the spectra is not well suited to study 
a modulation on this time-scale. Nevertheless, phasing the 
optical spectra on 1.23799 d does not provide a convincing 
result. 

The stellar wind changes are extensive over short-time 
scales (cf. Fig. 2). The Hell 4686 emission equivalent width 
changes by ~ 30% over 24 hours, and by ^ 10% in just 30 
minutes. Even in the latter rapid fluctuation the outflow is 
simultaneously altered over almost 300 km s~^ blueward and 
redward of line centre, i.e. ~ ± 0.25 of the wind terminal 
velocity in NGC 6826. Wc suggest that the stochastic nature 
of the dynamic spectra in Fig. 2 and the unstable frequencies 
evident in Fig. 4 point to the onset of clumping in the fast 
wind of the central star. 

However, Fourier analysis of the optical spectroscopy 
suggests at least one potentially interesting modulation time 
scale, as demonstrated in Fig. 7. We determined the central 
velocities of the CIV 580lA and 5812A emission lines by 
least squares Gaussian fits. The internal error bars on the 
velocities are < 1 km/s, and are "normalised" to the central 
velocity measurements of the Hp nebular emission. The am- 
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Figure 6. Time series of phase diagrams, constructed from some 
three years of LC Kepler photometry, folded around the subhar- 
monic frequency. Each phase diagram comprises a ~ 2-month 
subset of data. Time runs from top to bottom. Fits to the phase 
plots are shown as well. 
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plitude spectrum shows the strongest peak (fspec) at 7.821 
, or Pspec ~ 3.1 hr, with S/N — 3.9, just below our 
significance limit. 

A phase diagram relative to this period (middle panel 
of Fig. 7) folds nicely as to be expected. In the uppermost 
panel of Fig. 7, the dynamic spectra of the He 11 5411 A ab- 
sorption feature have been folded on the 3.1-hr period. The 
phase coverage is reasonable, and the image suggests some 
evidence for organised behaviour, such as the red-to-blue 
enhancement of the flux with respect to the mean. Unfortu- 
nately, this 3.1-hr spectral line modulation has no convinc- 
ing counterpart in the Kepler photometry, even when only 
analysing the parts of the light curve coinciding with the 
spectroscopic observations (Fig. 8). Although we cannot rule 
out a binary system based solely on the optical spectroscopy, 
we nevertheless do not detect a consistent periodicity in the 
radial velocities of different emission and absorption lines in 
NGC 6826. The spectroscopic time series does not fold well 
with a 0.619 or 1.238 d period, which is at least partly due 
to the sampling of the data. For the same reason a 3.1-hr 
signal cannot be seen in Fig. 2. 

3.3 Summary of observational results 

The variability of NGC 6826 is characterized by: 

• A basic 1.23799 ± 0.00001-d photometric modulation 
with a (mostly) double-wave variable hght curve shape and 
total light range not exceeding 0.01 magj. 

• Irregular variability with a light range of about 
0.05 mag'^ and a time scale of a few hours. 

• Optical stellar absorption and emission lines variable 
on time scales of a few hours, a possible periodicity of 3.1 
hr, with evidence for organized behaviour. 

• Discrete Absorption Components in wind-sensitive ul- 
traviolet spectral lines, with at least two episodes occurring 
over ~8.5hr (Prinja et al. 2012). 

• Red noise in the amplitude spectrum as a consequence 
of the irregular light variability. 

There is no evidence for a spectroscopic counterpart of the 
photometric 1.23799-d modulation, and no photometric ev- 
idence for the DAC re-occurrence time. Also, there is no 
correspondence between the optical radial velocity changes 
and the simultaneous Kepler photometry. 



4 DISCUSSION 

The lack of radial velocity variations with the same pe- 
riod as the only coherent signal in our data, its double- 
wave shape, and the temporal changes of this shape strongly 
argue against a pulsational or binary origin. The observed 
temporal behaviour of this signal is rather consistent with 
rotational modulation, i.e. we see areas of different surface 
brightness moving in and out of the line of sight. Prinja et 
al. (2012) measured wsini = 50 ± 10 kms"^ for NGC 6826. 

^ As estimated in Sect. 2.1, the intrinsic light amplitude would 
be a factor of four or more higher due to the presence of the 
nebula. 
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Figure 7. Bottom panel: amplitude spectrum of radial velocity 
shifts measured in the C IV 580lA and 5812A emission lines. 
Middle panel: the C IV central velocities phased on the 3.1-hr 
period and folded over two cycles. Top panel, upper half: the 
combined CFHT and NOT time series of the He II 5411 A ab- 
sorption feature shown as a dynamic spectrum now phased on 
the 3.1-hr period. Top panel, lower half: the mean He II 541lA 
line profile. 
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Figure 8. Simultaneous C IV radial velocities (upper panels) and Kepler photometry (lower panels) of NGC 6826. The non-overlapping 
parts of the light curves have been removed for easier comparison. The light and radial velocity variations do not correspond. 



With the temperature (Tes = 46 000 K) and surface grav- 
ity (logg = 3.8) of the central star spcctroscopically deter- 
mined by Kudritzki, Urbaneja & Puis (2006), these authors 
obtained R = 1.8 Rq. This gives a maximum rotation pe- 
riod of about 1.8 d. The period of 1.23799 d inferred above is 
consistent with this value, and implies an inclination angle 
of the rotational axis of the order of 45° . This leaves suffi- 
cient horizon for brightness inhomogeneities to appear and 
disappear. The shape of this photometric modulation and 
its slow changes may be interpreted as two "spots" whose 
size or surface brightness contrast somewhat vary with time. 

Concerning the hourly photometric variability, the time 
scales themselves immediately rule out a binary or rotational 
origin; a potential secondary component would be located 
within the primary CSPN, or the CSPN would needs to ro- 
tate faster than breakup speed, respectively. The absence of 
any short-term periodicity in the whole data set and the ir- 
regularity of these variations argue that we do not see a pho- 
tospheric phenomenon. This is also am argument that pulsar 
tions axe not observed in the Kepler light curves: pulsation 
occurs at distinct frequencies. Even if it was stochastic, it 
can only appear and disappear at stellar eigenmode frequen- 
cies. Therefore one would expect to see at least some peaks 
standing out in the periodogram of the SC data (Figs. 3), 
which is not the case. However, we cannot claim that pulsa- 
tions are not present at all because the photometry contains 
little flux that originates in the stellar photosphere. 

We evaluate the time scales of possible central star pul- 
sations. Prom the most suitable sequence of envelope models 
by Gautschy (1993) with respect to NGC 6826 we estimate 
pulsation constant Q = P ^ p/ PQ = 0.042 d for the funda- 
mental radial mode, that is excited by the K-mechanism in 
almost all of these models. This pulsation constant corre- 
sponds to a period of the radial fundamental mode of NGC 
6826 of 2.7 hr. Given tiie uucertaiutics in the stellar param- 
eters and models, this can only be seen as a rough estimate. 
We note that envelope models do not allow the computation 
of nonradial modes, but the presence of strange modes with 
similar periods as the radial pulsations was noticed in the 
models by Gautschy (1993). 

It is important to constrain under which circumstances 
we would have missed the detection of binary-induced vari- 
ability. A binary would reveal itself through radial velocity 



or light variability. The photometric noise level in our data 
is 0.072 mmag at worst. This translates into a detection 
threshold for periodic signals of 0.29 mmag. Given the di- 
lution of any given variability signal due to the presence of 
the nebula (Sect. 2.1), we conservatively estimate that we 
would have detected any periodic signal with an intrinsic 
amplitude exceeding 1.5 mmag in our data. 

Morris (1985) derived expressions for the amplitude of 
cUipsodial light variations. We applied these in combina- 
tion with the limb darkening and gravity brightening coeffi- 
cients for the Kepler bandpass by Claret & Bloemen (2011), 
and adopting a central star mass of 0.74 M© (Kudritzki 
et al. 2006). This resulted in the constraint that we could 
have detected binaries with orbital periods shorter than 
seven days under favourable inclinations. However, a stel- 
lar (M ^ 0.08 Mq) companion in a 1.238-d orbit would only 
have escaped our detection in combination with an orbital 
inclination below 26°. 

Under the assumption of a random orientation of the 
orbital plane in space, the probability that it is observed 
below a certain inclination angle is 



p{< i) = 1 — cosz. 



(1) 



Further assuming that the CSPN rotation axis is normal to 
this orbital place, this probability can be written as 



P(< i) = 



cos i) 



i ^ sin ^{{vaini)/vb) 



^vl-(vsini)^ (2) 

p{< i) = ^ i < sin~'^((t;sini)/t;6), 

where Vb is the stellar rotational break-up velocity (see 
Bernacca (1970) for a detailed discussion). 

We can estimate the breakup velocity of the CSPN by 
using the spectroscopic parameters quoted earlier and 



2ttR. 



3/2 



(3) 



where Req is the equatorial radius of the star. Following Reid 
et al. (1993), Req = 1.5 i?, where R is tiie stellar radius in 
the non-rotating case. This calculation yields Pcrit = 7.8 h, 
hence Vb = 280kms~^, hence sm~ '^{{v sin i)/vb) = 10.3°. 
With these values, probabilities calculated with Eq. 2 are 
only 1.6% higher as if they were derived with Eq. 1. Re- 
turning to the hypothesized stellar companion in an orbit 
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inclined by less than 26°discussed earlier, the probability 
that we missed its detection therefore becomes 10%. 

A stronger constraint arises from the lack of a detection 
of periodic radial velocity changes. We again used a central 
stax mass of 0.74 M©, and a generous upper limit for its 
periodic radial velocity variations of 3.3 kms~^ . For orbital 
periods shorter than seven days (the length of our spectro- 
scopic time scries), wc should have detected all stellar com- 
panions with an orbital inclination above 44°; for hypothe- 
sized companions exceeding 0.2 Mq, the limit on the orbital 
inclination is i < 31°, corresponding to a 14% probability of 
missing a detection. A stellar companion in a 1.238-d orbit 
would not have been detected in radial velocity with an or- 
bital inclination below 10° . However, in that case the CSPN 
rotation axis cannot be normal to the orbital plane because 
the central star would have to rotate above break-up speed. 

The referee remarked that the shape of the NGC 6826 
nebula and the presence of Fast-Low-Ionization Emission 
Regions (FLIERs) in it (BaUck et al. 1994) suggest that the 
orbital plane of a hypothetical binary causing these features 
would be at low inclination. Whereas such an interpretation 
cannot be ruled out completely, it needs to be reconciled 
with the orientation of the CSPN's rotation axis keeping in 
mind its critical rotation rate. Furthermore, such a low in- 
clination is inconsistent with our interpretation of the 1.238- 
d variation in terms of a rotational modulation originating 
from the central star, unless one postulates a rotation axis 
significantly deviating from the normal to the orbital plane. 



5 INTERPRETATION 

Following their analysis of one year of Kepler LC data on 
NGC 6826, Jevtic ct al. (2012) suggested that a combination 
of stellar pulsation and interaction with a close companion 
is a possible explanation of the central star's variability. Ac- 
cording to our results, such a scenario is very unlikely. Our 
measurements rather support another possibility mentioned 
by Jevtic et al. (2012), variable features on the central star's 
surface in combination with inhomogeneities and variations 
in the density structure of its wind. Such variability is obser- 
vationally well established in massive hot stars (e.g.. Fuller- 
ton 2003). 

Large-scale variations in the wind structure arc believed 
to be seated at the base of the wind, in so-called Co- Rotating 
Interaction Regions (CIRs, Mullan 1984). In brief, CIRs 
arise due to intensity variations on or near the stellar photo- 
sphere, such as spots or nonradial pulsation patterns. Such 
intensity or temperature variations modify the outflow ve- 
locity of the spherical wind. Hence, the parts of the wind 
originating at these locations are differently affected by the 
driving force, and collide with particles emitted from other 
surface regions. Due to stellar rotation spiral-shaped density 
fluctuations in the wind arise and cause the spectroscopically 
observed DACs. These DACs occur on time scales related 
to the stellar rotation period (several days) , and drift slowly 
through the line profiles with respect to the wind velocity. 
As exemplary literature wc refer to FuUerton et al. (1997) for 
an observational study and Cranmer & Owocki (1996) and 
Lobel & Blomme (2008) for hydrodynamical model compu- 
tations. 

In addition to the DACs, "modulations" have also been 



observed to move through the line profiles. They travel at 
a much faster rate than the DACs and have considerably 
shorter recurrence times (FuUerton et al. 1997). These "mod- 
ulations" alter the mean flux, as opposed to the DACs that 
are pure absorption features. Radiative transfer modelling 
(Lobel, Toala & Blomme 2011) implies that these axe only 
slightly bowed large-scale density enhancements and local 
wind velocity variations that radially protrude into the equa- 
torial wind. The density enhancement is of the order of 10%, 
and can result from mechanical wave action at the base of 
the stellar wind, such as produced by nonradial pulsations 
(e.g., see Kaufcr ct al. 2006). The photospheric cause of these 
variations is called Rotational Modulation Regions (RMRs). 

The variability of some Wolf-Rayet (WR) stars may be 
related to such variations, although they have been inter- 
preted slightly differently. Periodic photometric variations 
with time scales of days have been attributed to CIRs (e.g., 
Chene et al. 2011), and irregular short-term variations with 
time scales of hours have been reported in addition. Their 
spectroscopic counterparts are believed to be moving fea- 
tures in optical profiles of wind-sensitive lines, observed not 
only in WR stars, but also in Of stars (e.g., Lepine & Mof- 
fat 2008) and are interpreted as clumps propagating in the 
wind. 

The temporal behaviour of NGC 6826 conforms to these 
scenarios. Taking the photometric 1.23799 d period to be due 
to rotation, the double-wave shape may be a manifestation 
of two co-rotating features. The recurrence time for DACs in 
the UV lines of NGC 6826 is not well constrained; Prinja et 
al. (2012) reported the detection of two sequential features 
over ~ 8 hours. It is interesting to note however that we 
never see more than two strong CIRs in OB staxs that have 
been monitored for one rotation period. 

The time scales of short term light variations are consis- 
tent with RMRs or clumps moving within the wind. As our 
spectroscopy implies that the variations in the wind are not 
local, they are most likely attributable to large-scale fea- 
tures. Therefore they may be associated with RMRs. The 
difficulty with this picture is that the spectroscopic data of 
massive stars, in which RMRs are manifested, imply that 
they occur quite regularly, whereas the short-term photo- 
metric variability is irregular in our case. Consequently, 
some periodically triggered mass loss must have lost mem- 
ory by the time it reaches the depth in the wind that the 
photometry samples. 

Albeit the present observations do give some insight 
into the nature of the variability of NGC 6826 (and other 
ZZ Lcp stars), more observational evidence is needed. Al- 
though observationally hard to obtain due to the faintness 
of the targets, spectroscopic evidence for wind clumping or 
eventual RMR presence in CSPN should be sought. On the 
other hand, photometric evidence for light variability on the 
RMR and CIR time scales for massive stars, focusing on 
non-WR stars, would also be useful. 

In this context it is interesting to note that Blomme et 
al. (2011) analysed CoRoT light curves of three massive O- 
typc stars and also found some apparently incoherent vari- 
ability, in all three targets. The light curves, amplitude spec- 
tra, and time-frequency analyses of these stars phenomeno- 
logically closely resemble what we report for NGC 6826. 
One may therefore speculate that this kind of variability 
is present in all hot stars (Blomme et al. (2011) and Han- 
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dler et al. (2012) discussed it in connection with massive OB 
stars), and that it is related to stellar wind variations. 



6 SUMMARY 

The observed photometric variability of the central star of 
NGC 6826 consists of a periodic modulation with a time 
scale of 1.23799 d and irregular light variations on time scales 
of a few hours. Optical spectra imply line variability on a 
similar time scale, with a possible period of 3.1 hr. Ultravi- 
olet spectra show DACs with a re-occurrence time of about 
8.5 hr. The only possible correspondence between the pho- 
tometric and spectroscopic variability may be the time scale 
of the hourly changes. 

The periodic photometric variation is best explained 
by rotational modulation. The 1.23799 d-period we derived 
is similar to asteroseismic rotation periods of single white 
dwarf stars (e.g., Winget et al. 1991). Time-series photom- 
etry of ZZ Lep stars may therefore reveal their rotation pe- 
riods, but large data sets are needed to detect these peri- 
odicities due to the dominating short-term variability, that 
is due to changes in the stellar mass loss. This behaviour is 
similar to what has already been observed in massive OB 
and Wolf-Rayet stars, suggesting that the same mechanism 
may be responsible for the variations in all hot-star winds. 
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